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Nematic liquid crystals are well modeled as a fluid of rigid rods. Starting from this model, we 
use a Poisson-bracket formalism to derive the equations governing the dynamics of nematic liquid 
crystals. We treat the spin angular momentum density arising from the rotation of constituent 
molecules about their centers of mass as an independent field and derive equations for it, the mass 
density, the momentum density, and the nematic director. Our equations reduce to the original 
Leslie-Ericksen equations, including the inertial director term that is neglected in the hydrodynamic 
limit, only when the moment of inertia for angular momentum parallel to the director vanishes and 
when a dissipative coefficient favoring locking of the angular frequencies of director rotation and 
spin angular momentum diverges. Our equations reduce to the equations of nematohydrodynamics 
in the hydrodynamic limit but with dissipative coefficients that depend on the coefficient that must 
diverge to produce the Leslie-Ericksen equations. 

PACS numbers: PACS numbers: 61.30.-v,47.50.+d 



I. INTRODUCTION 

The Leslie-Ericksen (LE) equations p] for the dynam- 
ics of nematic liquid crystals have been a bulwark of liq- 
uid crystal science since they were first derived over forty 
years ago. For a historical account of their derivation see 
Ref. 0- They indisputably provide a correct theoretical 
description of the almost limitless variety of dynamical 
phenomena that nematic liquid crystals can exhibit, from 
simple shear flow to hydrodynamic instabilities to com- 
plex switching in display cells. 

The equations originally derived by Leslie and Erick- 
sen are not completely hydrodynamical. They contain 
an inertial term in the equation for the director, n, spec- 
ifying the direction of molecular alignment, that leads 
to modes that decay in microscopic times. When this 
term is ignored, the resulting equations are purely hy- 
drodynamical with mode frequencies that all vanish with 
vanishing wavenumber. Subsequent treatments jj, Q of 
the dynamics of nematics produced purely hydrodynam- 
ical equations from the outset. To our knowledge, all of 
the many experimental verifications of the validity of the 
LE equations probe only the hydrodynamic limit; they 
do not test the existence of or the form of the nonhydro- 
dynamic part of the original LE equations. 

The LE equations and purely hydrodynamic treat- 
ments of nematodynamics differ most profoundly in their 
treatment of kinetic energy. In the hydrodynamic theo- 
ries 0, Q , the momentum density g and its related ve- 
locity field v measure the momentum of all mass points 
in the medium including those along the full length of 
rigid mesogens. This momentum is a conserved quan- 
tity and is necessarily hydrodynamic. The energy den- 
sity is g 2 /(2p), where p is the mass density. Since g is 
the total momentum density, it contains all information 
about angular momentum, and it is not necessary to in- 
troduce additional variables to describe what we will call 



the spin angular momentum associated with rotation of 
constituent rigid molecules about their centers of mass 
0, H • The other variables that appear in the hydrody- 
namic treatment are the conserved mass density and the 
Frank director n. In the LE treatment, there are two 
contributions to the kinetic energy: a translational part 
g 2 /(2p) and a rotational part, I±(n x h) 2 /2, where I± 
is a moment of inertia density, arising from motion of 
the director. The interpretation of this decomposition 
of the kinetic energy is that g is now the center-of- 
mass momentum density, which is a conserved variable, 
and that the director contribution to the kinetic energy 
arises from the spin angular momentum. There are now 
two contributions to the angular momentum: the spin 
angular momentum and that arising from center-of-mass 
motion. Neither contribution is individually conserved, 
but their sum is. In the LE treatment, the equation for 
the director is basically an equation for the spin angular 
momentum, which is neither a conserved nor a hydrody- 
namical variable and which, therefore, has nonhydrody- 
namic decays in it. 

Though the LE equations are internally consistent and 
reduce to the correct hydrodynamical form when the 
inertial term is ignored, they in fact do not provide a 
correct description of spin angular momentum. Spin 
angular momentum is an independent dynamical vari- 
able that is not locked to the director, though in steady 
state situations it does relax to a value determined by 
the local director and its rotation rate. The spin ki- 
netic energy is determined by the spin angular momen- 
tum and not by the dynamics of the director. In this 
paper, we describe the dynamics of nematics in terms 
of their non-conserved spin angular momentum density, 
their conserved mass density and center-of-mass momen- 
tum density, and their director. For simplicity, we con- 
sider isothermal processes only, and we ignore the equa- 
tions of energy conservation. We use the Poisson-bracket 
approach [E S 0, S E3, 113, E3, E3, EE| to derive the equa- 
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tions of motion for these variables. The hydrodynamical 
limit of our equations is identical to that of the LE equa- 
tions but with a slightly different interpretation of some 
dissipative coefficients. Our equations also reduce to the 
full LE equations when an appropriate viscosity diverges 
and spin angular momentum parallel to the director is 
ignored. Previous investigations [l3l H(| have noted that 
spin angular momentum should be treated as an inde- 
pendent variable and argued that it will decay in micro- 
scopic times to a value determined by the director. They 
do not, however, provide a detailed prescription for how 
this decay occurs or the conditions under which the orig- 
inal LE equation can be retrieved. Finally, our equations 
reduce to the hydrodynamical equations for rigid rotors 
on a lattice pH | when coupling to center-of-mass motion 
is turned off. 

Since the derivation of our results is at times alge- 
braically tedious, we review our main results in Sec. ITTI 
We present first the dynamical equations that result from 
our analysis and show how they reduce to the LE equa- 
tions and to true hydrodynamical equations in the appro- 
priate limit. In Sec. lIIII we briefly summarize the Poisson- 
bracket formalism. In Sec. II VI we introduce the fluid of 
rigid rods and the relevant dynamic variables and their 
Poisson brackets. Finally, in Sec. |Vl we derive the equa- 
tions of nematodynamics with spin angular momentum. 



II. REVIEW OF RESULTS 

We model our nematic as a collection of uniaxial rigid 
rods. The coarse-grained variables describing this system 
are the mass density p, the center-of-mass momentum 
density g — pv, the spin angular momentum density I, 
and the nematic director n. The angular momentum 
density I = ISl can be expressed in terms of a moment- 
of-inertia density = Tumrij + X±(5ij — nirij) and an 
angular frequency ft. The full equations for all of these 
variables are 
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h is the molecular field with com- 
nirij)5F/Sni, where F is the free 
energy, p is the pressure, <7y is the dissipative and direc- 
tor part of the stress tensor, A is the the symmetrized 



where AI = I\\ — I± 
ponents hi — — 



strain rate tensor with components (diVj + djVi)/2, and 
u> = (V x v)/2 is half the local vorticity. In these equa- 
tions, dA/dt — dA/dt + v ■ 'VA is the total derivative of 
any field A, cy = n ■ c is the component of any vector 
c along n and c± = n x (c x n) its component perpen- 
dicular to n, and An is the product of a matrix with 
a vector with components AijUij. The stress tensor cr^ 
can be decomposed into elastic, viscous symmetric, and 
viscous anti-symmetric parts: 



ij = ij + °i] + a ij i 



where 

e _ df 



a 



i j 



(2.6) 



(2.7) 



a ij = a i]kiA k i + -T A (e l ikn ] n l + e j i k n i ni){Q. k - ^k), 

(2.s; 



4' = 2 e ^ r ^-^) 
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with T^i = T^niTij + T^_(Sij — n.irij) and a fourth- 
rank tensor of uniaxial symmetry given explicitly in Eq. 

Equations (|2.1|) to (|2.9|) provide a complete description 
of the dynamics of nematics. Equations l|2.1() and 12.21) 
are the familiar conservation laws for mass and momen- 
tum. Equation (|2.3|l is the equation of motion for the 
director. It is similar to that of the full hydrodynamical 
theory derived by the Harvard group Q except that the 
reactive term f2 x n depends only on the spin frequency 
n and not on A and uj. If 7 -1 = 0, the director sim- 
ply rotates like a rigid-body axis perpendicular to the 
angular spin f2. When 7 -1 is non zero, director motion 
decays to f2 x n in a time of order K/^q 2 ) where K is 
a Frank clastic constant and q is the spatial wavenum- 
ber of director distortions. Equations (|2.4D and (|2.5|l are 
the torque equations for spin angular momentum. The 
spin frequency parallel to the director, described by 
Eq. H2.5f) . is a nonhydrodynamic variable that decays in 
microscopic times to wy plus nonlinear terms. It cannot 
be ignored except in the hydrodynamic limit or in the 
limit in which only the component of spin angular mo- 
mentum perpendicular to the director survives, i.e., when 
X|| =0, which occurs in the limit of perfect nematic order 
in a system composed of infinitely thin rods. Equation 
(|2.5|) contains a nonlinear term proportional to f2||f2j_ 
that has physical significance. As reviewed in Appendix 
1X1 it causes the axis with the highest moment of inertia 
to align along the direction of the angular momentum 
in rigid body motion with conserved angular momentum 
but not energy. The right hand sides of Eqs. I|2.4|l and 
(|2.5() are the time rate of change, dl/dt, of the spin angu- 
lar momentum projected, respectively, along directions 
perpendicular and parallel to n. In the absence of dissi- 
pation, dl/dt is simply the torque density nxh appearing 
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in Eq. I(2.4|) . Rigid-body rotation in which f2 = ui and 
A = is a stationary state in which dl/dt — 0. HQ ^ us, 
dissipative torques, given by the P^ and terms in Eq. 
(|2.4|) and ((2.5(1 . drive towards u;. Spin angular momen- 
tum is also reoriented via the T A term in Eq. QTty by the 
symmetric strain rate A when it is nonzero. The stress 
tensor contains a couple of terms not found in isotropic 
fluids. The anti-symmetric parts of the stress tensor pro- 
portional to r^4 and are dictated by the requirement 
that the total spin and center-of-mass angular momen- 
tum is conserved (see Appendix[5|) . The T A contribution 
to the symmetric part of the dissipative srress tensor afj 
is a consequence of an Onsager relation. 

When rjp T^, and T A are zero, spin angular mo- 
mentum is conserved, and additional diffusive dissipative 
terms proportional to V 2 Q must be added to Eqs. (|2.4|) 
and H2.5|) for a complete description. In this limit, Eq. 
(12. 3|) along with Eqs. (|2.4|l and (|2.5|) provide a hydrody- 
namical description of rigid rotors on a rigid lattice with 
frictionless bearings, which exhibit spin- wave excitations 
with a frequency dispersion lo ~ q with q the wavenumbcr 

M 

We can now consider under what conditions our equa- 
tions reduce to the original LE equations and how the 
hydrodynamical limit is obtained. We begin with ob- 
taining the LE equations. To approach the LE limit, we 
use Eq. ((2.3(1 to replace Qxn by (dn/dt) — r y~ 1 h T , where 
h T = h + AJOiifix- This converts Eq. to 



2±n x 



<ll ! 

d 2 n 1 dh 1 



dt 2 7 dt 
= an x (h — 7iAT — "f 2 An) 



In 0,1 — 

11 11 dt 



dn p?at 



7 



where we introduced 
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dt 



and 
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-flwn x q 
(2.10) 

(2.11) 



^n7i = -A71 



r 



(2.12) 

Equation (12.10(1 reduces identically to the original LE 
equation lj for the director with left-hand side equal to 
X\n x d 2 n/dt 2 , a = 1 and 71 = P^ 2 , when 7 — > 00 and 
X|| = 0. The first condition, 7 — > 00 constrains dn/dt 
to be r2 x n. The second condition, Zj| = 0, is equiv- 
alent to there being no rotational kinetic energy associ- 
ated with 0|| and is one that is tacitly assumed in the 
original LE approach for which the spin kinetic energy 
density is I±(n x n) 2 /2. We will show that T|| vanishes 
for rigid rods when the Maier-Saupe order parameter S 
equals one, i.e., only when there is perfect order. If ne- 
matogens are modelled by more complex rigid structures 
than thin rods, Iu would be nonzero even for S = 1. 
In the LE limit, Eq. ((2.5(1 for f2|| implies that f2ii = uu. 
Together with Q± — uj± = n x TV, which follows from 
dn/dt = ft x n and the definition for AT, the equations 



for the stress tensor [Eqs. 1(2.8(1 and 1(2.9(1 ] assume exactly 
the form of the LE stress tensor. Thus to reiterate, the 
LE equations describe a nematic liquid crystal in which 
the director is forced to follow J~2 x n and the moment of 
inertial density parallel to the director is zero. Neither 
of these conditions apply in general. 

To obtain the hydrodynamic limit, we discard all terms 
that are higher order in time and space derivatives than 
the dominant ones. This means that we can ignore the 
dVL\\/dt and the nonlinear terms in Eq. ((2.5|) relative 
to f2|| — Thus to hydrodynamic order, we can set 
fin = oj\\. This procedure effectively removes fin from 
the problem. Similarly, we can ignore the 0|| term in the 
director equation [Eq. 1(2.3(1 ] and all of the terms on the 
left-hand side of Eq. ((2.10(1 . The latter condition gives 
the familiar LE equation, h — 71 N + 72 (An)±, for the 
director in which the inertial term is neglected. To ob- 
tain the hydrodynamic limit for the stress tensor, we use 
the hydrodynamic limit of Eqs. 1(2.3(1 and 1(2.10(1 and the 
relations in Eq. 12.12(1 to set 



ill u>> x n = -^{An)± Cl.YX) 



and obtain 



vt! = 7;li{ n i N j ~ n i N i) + ^72[ni(An)±j - nj(An)±i] 



1 



(2.14) 



a ij = otijkiAki + -^l2{riiNj + rijNi) 



A\2 



1 {T A ) 



-[m(An)xj +rij(An)xi]. (2.15) 



2 7 + P 

This is precisely the LE stress tensor in the hydrody- 
namic limit. 

An important consequence of this analysis that treats 
angular momentum as an independent variable is that it 
demonstrates that two distinct effects contribute to the 
viscosity 71: the director damping measured by 7 and 
the rotational friction measured by . 71 is the parallel 
combination of 7 and [see Eq. ((2.121) ]. 

In Ref. ^3 we derived dynamical equations for the full 
nematic order parameter Q, also called alignment tensor, 
using the Poisson-bracket formalism without, however, 
introducing the spin angular momentum density as a sep- 
arate dynamic variable. With the approach presented in 
the following, we could also derive dynamic equations 
for Q and then by projection on the uniaxial part of Q 
arrive at an additional dynamic equation for the scalar 
order parameter. A similar consideration following Eqs. 
(12. 101) - ((2. 12(1 should then lead to the extended LE equa- 
tions of Ref. [l2j, where a variable S is taken into account. 
The field S is a nonhydrodynamic variable that relaxes 
in microscopic times. It does not, however, contribute to 
dissipative coefficients in the hydrodynamic limit as the 
spin angular momentum, which also has a rapidly decay- 
ing nonhydrodynamic component, does [see Eq. 1(2.12(1 ]. 
We will, therefore, not treat S (or the biaxial part of Q) 
in what follows. 
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In the remainder of the article we give a detailed ac- 
count of how the set of equations discussed in this sec- 
tion were derived. We start with a short review of the 
Poisson-bracket formalism. 



III. GENERAL FORMALISM 

In this section we collect the important formulas of 
the Poisson-bracket formalism. A more thorough expla- 
nation including original references can be found in our 
previous article 0] and in |T(1 |. 

We consider a systems whose microscopic dynam- 
ics is determined by canonically conjugate variables q a 
and ir a for each particle a and a microscopic Hamil- 
tonian H({q a }, {n a }). Rotational degrees of freedom 
may be included in the coordinates q a and momenta 
iz a . We are interested in the slow dynamics of of a 
set of macroscopic field variables ^(a:, t) (// = 1,2,...) 

obtained from microscopic fields ^^(x, {q a }, {tt 01 }) by 
coarse-graining over spatial fluctuations on the micro- 
scopic level; $ M (a;,t) = [^^(x, {q a }, {ir°}) ] c , where the 
symbol [. . .] c specifies the coarse-grained averages. The 
statistical mechanics of the macroscopic fields t) is 

determined by the coarse-grained Hamiltonian 7i[{$ M }]. 

Following the theory of kinetic or stochastic equations, 
the macroscopic variables evolve according to 



sn 



at ' ' ~ v ^ x) r ^ 5$ v {x) ' 



(3.1) 



where we disregard any noise. The reactive term V^(x), 
also called the non-dissipative or streaming velocity, is 
expressed with the help of Poisson brackets as 



V^x) = - J d 3 x'V„ u (x,x 



,. sn 



(3.2) 



where Einstein's convention on repeated indices is under- 
stood and 

P iw {x,x t ) = {$»{x),$ v {x')} = -V^{x',x) (3.3) 

denotes the Poisson bracket of the coarse-grained vari- 
ables. It is defined as the coarse-grained average of the 
microscopic Poisson bracket: 

*„(*')} = [{$„(*), $„(*') }]c , (3.4) 

where flTl 



force 6H/8$> v (x), which together with &v(x) forms a pair 
of conjugate variables. The dissipative tensor r M „ may 
depend on the fields $ M and it may also contain terms 
proportional to —V 2 . It is determined by three princi- 
ples. First, the dissipative contributions to the equation 
for dQ^/dt must have the same sign under time reversal 
as 3"^ (and thus the opposite sign to that of dQ^/dt). 
Second, T ^ has to reflect the local point group symm- 
metry of the dynamical system, and third, it has to be a 
symmetric tensor at zero magnetic field to obey the On- 
sager principle |l8j |. In the following, the last point will 
be important in identifying the proper dissipative terms 
in the momentum balance. 



IV. POISSON BRACKETS FOR NEMATIC 
LIQUID CRYSTALS 

A. Model molecule and dynamic variables 

We model our system as a fluid of uniaxial rigid rods 
of length a and mass m. We describe the position of 
molecule a by its center-of-mass coordinate x a and its 
orientation by the unit vector u a . The center-of-mass 
momentum is p a = mv a = mx a where dot means total 
time derivative. With the help of the molecular order- 
parameter tensor Q a with components 



Qi 



-8 ■ 



(4.1) 



the molecular moment-of-incrtia tensor (relative to the 
center of mass) reads 

Jg = I^vfvf + I^ - vfvf) = AJQg + (§/l + i/y)% 

(4.2) 

where I± and 7m < I± are the moments of inertia for 
rotations about axes perpendicular and parallel to i> a 1 
respectively, and AI = Li — I±. Note that Iu = 0, and 
AI = —I± in the limit of infinitely thin rods. The tensor 
I a relates the spin angular momentum of a molecule, l a , 
to its angular velocity fJ Q : 



I? 



i a n a 

n 3 



(4.3) 



Note that l a is always perpendicular to i> a as it should 
be for an infinitly thin rod. 

The microscopic Poisson bracket in Eq. (|3.5|l consists 
of a spatial and an angular part. The spatial contribution 
arises from the coordinate x a and its conjugate momen- 
tum p a , which fulfill the canonical Poisson bracket [Tt| : 



E 



d$ M (x) d$ u (x') d^^x) d$ v {x') 



dqf 



(3.5) 



Since we only employ a restricted number of macro- 
scopic variables, all the "neglected" microscopic degrees 
of freedom give rise to the dissipative term in the kinetic 
equation (|3.1|) that is proportional to the generalized 



\p?,x?} = 8 a P8 i: 



(4.4) 



where 8 a ^ and 5y are Kronecker symbols. However, the 
unit vector u a and the angular momentum l a are not 
canonically conjugate to each other since their Poisson 
bracket pj| 



(4.5) 
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does not have the canonical form (sijk denotes the Levi- 
Civita symbol). We could now introduce appropriate 
pairs of conjugate angular coordinates and momenta via 
a microscopic Legendre function. Instead, we follow an 
alternative route. It turns out that, in the following, Eq. 
(|4.5() and the additional formula 



{if,#} = -* tt/ W£ 



(4.6) 



are sufficient for calculating the angular part of Poisson 
bracket ((3.5|) . All other Poisson brackets, in particular 
the ones between the spatial and angular variables, are 
zero. The Poisson brackets in Eqs. 1(4. 5fl and l|4. (jfl are 
a consequence of the fact that If is the generator of ro- 
tations about the molecular center of mass 0. They 
can, or course, be derived from a Hamiltonian formalism 
for rigid rods in which the rotational kinetic energy de- 
pends on two Euler angles for the case of infinitely thin 
rods (I\ | = 0) and three Euler angles for the more general 
case. 

We are now ready to define the relevant microscopic 
field variables and their coarse-grained counterparts. The 
conventional microscopic definition of the density of mass 
and center-of-mass momentum are: 



p(x) — m S(x — x a ) 

a 

g(x) = 5><M(*-< 



(4.7) 
(4.8) 



which, after coarse graining, result in the macroscopic 
variables p(x) = \p(x)] c and g{x) = [g(a;)] c = p(x)v(x). 
The last term defines the macroscopic velocity field v(x). 
In a similar manner we introduce the macroscopic field 
Q(x) of the nematic tensorial order parameter pL Hg . ll9l 
|2C[ . also called the alignment tensor [2l| : 



p(x) 



Q(x) = 



^Q a 8{x-x a ) 



(4.9) 



using Q a from Eq. I|4.1|l . The factor p{x)/m is intro- 
duced to make Q(x) unitless. With the microscopic def- 
inition for the density of the moment-of-inertia tensor, 

a. 

= £[AIQg + (§/_l + §/||)<y<5(* - * a ) , 

(4.10) 

we obtain the coarse-grained moment-of-inertia density 
p{x) 



Tij (x) 



(4.11) 



Finally, the microscopic field of the density of spin angu- 
lar momentum is 



T(x) = ^l a 5{x~x a ) 



(4.12) 



Its associated coarse-grained variable is l{x) = [l(x)] c = 
T(x)n(x), where the last expression defines the macro- 
scopic field tt(x) of angular velocity in full analogy to 
v(x). 

We are interested in the dynamics of the nematic phase 
where the orientational order is uniaxial. The alignment 
tensor therefore assumes the form 



i(x) = S[ni{x)nj(x) 



(4.13) 



where, on average, the molecules point along the director 
n(x). The Maier-Saupe order parameter 5* is constant 
in the nematic phase. With the uniaxial Q{x) of Eq. 
(14.13(1 . the moment-of-inertia density T(x) of Eq. ((4.11(1 
becomes 



(4.14) 



with 



1\ 



m \3 
p (2 



m V 3 



Iais) 



Its anisotropy is quantified by 



A2 = lu -1_l = AI—S 
m 



(4.15) 



(4.16) 



Note, as indicated, in Sec. |H] that In = in the limit of 
infinitely thin rods (7|| — 0) and perfect nematic order 
(5 = 1). 

Thus, our set of dynamic variables is {p, n,g,l} for 
which we have to determine all possible Poisson brackets. 



B. Poisson Brackets 

The calculation of the Poisson brackets is straightfor- 
ward. In addition to the comments about the angular 
variables in the previous section and the antisymmetry 
relation expressed in Eq. ((3.3(1 . we use properties of the 
5 function summarized as 



S(x - x') 
f(x)S(x-x') 
V. t 5(x - x 1 ) 



5{x' - x) 
f(x')S(x-x') 
-%5{x x') 



(4.17a) 
(4.17b) 
(4.17c) 



where = d/dxi, V- = d/dx^, and f(x) is an arbitrary 
function including the 5 function itself. 

In the following, we list all the non-zero Poisson brack- 
ets which determine the non-dissipative velocities of our 
dynamic variables. The dynamics of the center-of-mass 
density p(x) is provided by 



{p{x),g i {x')} = V i 5(x-x')p{x') 



(4.18) 



To derive the Poisson brackets of the director, we first 
calculate the Poisson brackets of the alignment tensor. 
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According to the definition Q4.9[l . we only have a micro- 
scopic expression for p(x)Q(x) but not for Q(x) alone. 
An analogous, however more c omp licated, situation oc- 
curred in our previous article To calculate, e.g., 
{Qij(x),gk(x')} : we apply the Product rule for Poisson 
brackets to {p(x)Qij(x), gk(x')} and arrive at 

1 



{Qij{x),g k {x')} = 



p(x 



■{p(x)Q ij (x),g k (x')} 



Qij 

p{x) 



{p(x),g k (x')}4.l9) 



The first term on the right-hand side and the second 
term, already known from Eq. I|4.18|) . then combine to 
yield 



{Qij(x),g k (x r )} = V k 5(x - x')Q ij (x') . 
In the same manner, we calculate 
{Qij(x),l k (x r )} = -\e ijk Q a {x) +eikiQji{x)]5(x 

where we have used the product rule and Eq. 
evaluate the microscopic Poisson bracket {Qfj,l k } and 
the fact that {p(x),l k (x')} = 0. The Poisson brackets for 
the director now follow by projection from the uniaxial 
representation (|4. of the alignment tensor (see Ref. 



(4.20) 



x ) , 
(4.21) 

m to 



repi 

m 



for details): 



{n i (x),g j (x')} = ^S&{Qu(x),g j (x')}n l (x) , (4.22) 



where 



(4.23) 



is the projector on the space perpendicular to n{x). The 
same formula is valid with gj(x') replaced by lj(x') so 
that Eqs. (|4~2UI) - (|4~23j) finally give 



{ni(x),gj(x')} 
{n t (x),lj(x')} 



[Vjni(x)]6(x - x') (4.24a) 
-e ljk n k {x)5(x - x') . (4.24b) 



Note that the Poisson brackets of Qij(x) and rii(x) with 
g k (x') are much simpler in this formulation with spin 
ang ular momentum than they are in the alternative one 
|llL |22j in which there is no spin angular momentum and 
g is the total rather than the center-of-mass momentum 
density. In particular, the director-momentum bracket 
\ij k V k 5(x — x 1 ) that plays such an important role in the 
latter formulation is not present in the current one. 

The reactive velocity of the translational momentum 
follows from 



{gi(x),p(x')} 
{g t (x), nj (x')} 
{9i(x),9j{x')} 



p{x)Vi8{x - x') (4.25a) 
-[V l n J {x)]8{x -x') (4.25b) 
-V{[5(* - x') 9j (x')] 
+Vj5(x - x')g t (x') (4.25c) 
{gtix),^')} = IjWViSix - x') . (4.25d) 



Eqs. (|4.25a|l and l|4.25bf) are related to Eqs. (|4.18(l and 
(|4.24a|l by the antisymmetry relation of the Poisson 
brackets whereas Eqs. (I4.25c|l and (|4.25d|) are readily 
calculated. Again the missing term in the momentum- 
director bracket compared to Ref. ^l| is compensated 
by the additional Poisson bracket il4.25d(l . 

Finally, the non-dissipative dynamics of the angular- 
momentum density is governed by 



{h(x),gj(x')} 
{h(x)Mx')} 



-e ijk n k (x)8(x - x') (4.26a) 
li{x')Vj5{x - x') (4.26b) 
-e ijk l k (x')6(x - x') . (4.26c) 



V. NEMATODYNAMICS WITH SPIN 
ANGULAR MOMENTUM 

Following the systematic structure of the theory out- 
lined in Sec. IIIII we now derive the full set of equations 
as presented and discussed in Sec. [H] We first calculate 
the reactive and dissipative velocities needed to formu- 
late the dynamic eqations for the set of dynamic variables 
{p, n,g, 1} and then introduce the spin angular velocity 

n. 



A. Non-dissipative velocities 

To calculate the non-dissipative velocities from Eq. 
(|3.2|) . we need the Hamiltonian 



n = 



^Im^Mx) 

+F[p(x),n(x)\ . 



d 3 x 



(5.1) 



It consists of a kinetic part, subdivided into a 
translational and rotational term, and a free energy 
F[p(x), n{x)\ = J f(p, n, Vn)d 3 i, which is Frank's free 
energy plus a term depending only on p. In the following, 
we will need derivatives of the inverse of the moment-of- 
inertia density such as dX~^ jdy where y stands for p or 
n. Taking the derivative of XifL~^ = Si k with respect to 
y, we find 



dy 



17/ = 



d 







(5.2) 



The non-dissipative velocity for the density is simply 
V p = —V • g(x) which immediately gives the mass- 
conservation law. For the director, we use the Poisson 
brackets (I4.24|) and the fact that STi/Sl^x 1 ) — 
rij(x') and arrive at 



1-^l k 



V? = -v(x) ■ Vm(x) + e ijk Vt j {x)n k {x) 



(5.3) 



The first term on the right-hand side is the convective 
derivative of n. The second term introduces a reactive 
coupling to the angular velocity Q(x). In the purely 
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hydroydnamical model 0, > this term is replaced by 
Xijk^ ' jVk coupling drii/dt to the symmetric and anti- 
symmetric parts of the deformation-rate tensor ViVj 

The non-disspative velocity V 9 is calculated with 
the help of the Poisson brackets H4.25J1 . Applied to 
the rotational part of the kinetic energy in Ti., i.e., 
^li(x)Xn (x)lj(x), they produce a contribution Vf ot to 
the non-dissipative velocity V 9 whose terms add up to 
zero. Specifically, we find 



V? 



-pVi(-fc 



+(V<n i )(-i fc Z l -g ; Jf-) -IjVify , (5.4) 



dp 

1 , , C^M 1 



which we rewrite as 



g _ 



fir- 1 / / 
! \^2 fci k 



~ Vi ( p 2 lkh ~dp L ') -ViVify) ■ ( 5 - 5 ) 

Using Eq. Q5.2J1 to evaluate the derivative with respect to 
p and the fact that Iki is linear in p [see Eqs. (|4.14ll and 
(|4.15() ]. one shows immediately that the three terms add 
up to zero, i. e., Vf ot — 0. All the other contributions to 
V 9 can be written in a compact form, 



V? = -V, 



9i{x)9j{x) 



p(x) 



- V t p + Vj<7. 



3 U V 



(5.6) 



as explained in Ref . • The first term on the right-hand 
side introduces the momentum flux tensor, the second 
term contains the pressure p, and cr E is the Erickscn 
stress tensor: 



Of 



(5.7) 



Note that compared to Ref. llj a term that contains the 
molecular field SF/5n is completely missing. 

Finally the Poisson brackets (|4.26(l give the reactive 
velocity of the angular-momentum density: 

rrl , sf 1 , , dI rs SF \ 

-Vj[k(x)vj(x)] + £ij k Slj(x)l k (x) . (5.8) 

Again, one can show using Eqs. lETHE)) and lt5~2f that the 
first and the fourth term on the right-hand side resulting 
from the rotational kinetic energy cancel each other so 
that we obtain 



V- = -Vj[h(x)vj(x)] +£ij k n k {x) 



SF 
Snj(x) 



(5.9) 



The first term introduces the angular momentum flux 
tensor in full analogy to the linear momentum and the 
second term is a coupling to the molecular field. 



B. Dissipative velocities and final equations 

1. Center- of -mass density 

For the conserved center-of-mass density, dissipative 
velocities are not allowed, and the mass-conservation law 
follows: 



2. Director 



(5.10) 



The time derivative dn/dt couples dissipatively only 
to forces conjugate to fields, n and p with the same sign 
under time reversal as n. A dissipative term proportional 
to n5Ti./dp, which has the correct sign under time rever- 
sal, cannot occur because it is always perpendicular to 
dn/dt 11]. A second dissipative term introduces a cou- 
pling to SH/Srii with a dissipative tensor <$^/7, where the 
projector Sf^ defined in Eq. I|4.23|l ensures that dn/dt is 
perpendicular to n and 7 is a rotational viscosity. To- 
gether with Eq. (|5.3() . the dynamic equation for the di- 
rector then reads 

drii „ n 1 T /l, , dlZ, 1 SF\ 

- W = -v Vn, + e ijk Cl jnk - -5 l3 {^hh^- + ^ J ■ 

(5.11) 

With the definition of the components of the molecular 
field, hi = —5fjSF/5rij, and the expression 



■ Xj_)S1||S1j 



(5.12) 



which we derive with the help of Eqs. 14.14fl and 15. 2|) 
and U — XijQj, we finally arrive at the director equation 
as presented in Eq. (|2.3(l in Sec. [H] 



3. Spin angular momentum density 

The fields SH/Sh = fi^ and 5H/Sgi = Vi have the 
same sign under time reversal as I and can contribute 
terms of the equation for dl/dt. To determine the form 
of these terms, it is important to realize that I and ft 
are pseudo vectors that do not change sign under space 
inversion but that the momentum density, g, and the 
velocity, v, are vectors that do change sign under space 
inversion. Thus a term directly proportional to Vi in the 
equation for dlt/dt is prohibited but one proportional to 
£!i is not. Pseudo vectors that are even under n — > — n 
can be constructed from the spatial derivatives of v and 
the director. Thus we look for dissipative terms con- 
taining these pseudovector combinations of Vi«j. The 
first pseudovector is u>i — e^feV jVk/2. Together with 
the angular velocity, it gives rise to the dissipative term 
— r^(Oj — uij) whose form is dictated by the require- 
ment that during a uniform rotation of the whole sample 
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(ft = lj) no energy is dissipated. A second pseudo vector 
of the velocity, which preserves the n —> —n symmetry 
of the nematic phase, is \(£ijinirik + Eikinirij)Ajk, where 
Ajk = (Vjffe + Vfci>j)/2 stands for the symmetrized ve- 
locity gradient. Furthermore, the third-rank tensor in 
front of Ajk is symmetric in j and k which is important 
for the next paragraph where we use the Onsager prin- 
ciple to find the dissipative velocities for the momentum 
density. Introducing the dissipative torque 



t[ = -Tfjflj - ujj) —(e lj inin k + £ ik ininj)A lk 



where 



r,, riiUj + Tj_(8ij - Turij) 



(5.13) 



(5.14) 



obeys the uniaxial symmetry of the nematic phase, and 
combining it with Eq. Ij5.9|l . we arrive at the formula 
describing the dynamics of I: 



— = -W j (l i v j ) + £ij k n k j^ + r t 



(5.15) 



To replace I by the angular velocity, we write 

I =J ll fl ll n+l ± fl± . (5.16) 

The time derivative of I involves d2 a /dt, where a repre- 
sents II or _L We find 



dig 

at 



dl a dp 
dp dt 



-Vi(X a Vi 



(5.17) 



where we used the fact that T a linearly depends on p 
[see Eq. Q4.15p ] and where we also employed the mass- 
conservation law dp/dt = —Vi(pvi). With Eq. (|5.17|) 
and the definition of the total time derivative, d/dt = 
d/dt + V{Vi, it is straightforward to show that 



dk dfl ± 



^ dn dfl\\ 



We introduce this term into the balance equation (|5.15(l 
for I and project it on n and the plane perpendicular to 
n to finally arrive at the respective Eqs. I|2.5|l and (|2.4|) 
in our review section [n] In deriving the last two terms 
of Eq. 12. 5|) . we also used n ■ dfl±/dt = —(dn/dt) ■ 
and replaced dn/dt by Eq. (|2.3|1 



4- Center-of-mass momentum density 

The dissipative term of the momentum balance is de- 
termined by the viscous stress tensor cr' which couples 
again to fli and ViVj as in the previous paragraph. The 
form of the dissipative part of the stress tensor is sub- 
ject to restrictions. First of all, because the total angular 
momentum (spin plus center-of-mass) is conserved, the 
antisymmetric part cr A ' of the viscous stress tensor and 



the dissipative torque r' of Eq. I|5.13|l are related (see 
Appendix [H| : 



1 



1 

2 



" 2' 

= o £ i3k^kli^l - Ml) 



+ -T A (njA in n n - mAj n n n ) . (5.19) 

To construct the symmetric part <t Si of the viscous stress 
tensor, we use Onsager's principle. It says that the dissi- 
pative fluxes r', cr s ' are coupled to the generalized forces 
ft — us, A by a symmetric, dissipative tensor. In symbolic 
notation this means 

-t' \ _( T n T A enn \ ( ft - u 
r 5 ' J ~ ^T^enn)' a ) \ A 

(5.20) 

where the superscript t in (enn) 1 stands for the appro- 
priately transposed third-rank tensor enn of Eq. I|5.13|l . 
The first line of the tensor equation reproduces the dis- 
sipative torque l|5.13|l . the second line gives 

1 



ctijkiA k i + -V (eukrijUi + eji k nini)(flk - uj k ) 



(5.21) 

where the viscosity tensor ot has the usual form required 
by the uniaxial symmetry of our medium (see, e. g., Ref. 
Ell): 

otijkt = ainiUjUkni + ^-(5 ik 5ji + 5u8 ]k ) 



a 5 + a 6 



(n^kSji + UjUkSu + riiniSjk + nj-nidik) 



(5.22) 



Adding up reactive and dissipative terms, the momentum 
balance finally reads 



dt 



-V j (^)-V i p + V J (4 + 4' + 4') . (5.23) 



The complete set of equations of nematodynamics in- 
cluding the spin angular momentum is reproduced in Eqs. 
(|2.1|) to H2.9fl . The Leslie-Ericksen and hydroydnamic 
limits of these equations are derived in Sec. 1TT1 
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APPENDIX A: NONLINEARITIES IN 
RIGID-BODY MOTION 

Our equation for the director reduces to that for the 
anisotropy axis n for a single uniaxial rigid body if n(x) 
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is independent of x and h = 0. If = and T A = 0, APPENDIX B: TORQUES AND STRESS TENSOR 

spin angular momentum is conserved [See Eq. (|5.15l) ]. 
i.e., dL/dt = d/dt(J d 3 xl) = 0, and Q\\ = L\\/I\\ and 
£l± = L±/I± where In ± = J d 3 xl\\^. The equation of 
motion for the anisotropy axis is then 



dn „ AI 
—— =S!xnH 

dt 7 



f2||fij 



(Al) 



This equation has a dissipative term implying that en- 
ergy is not conserved even though spin angular mo- 
mentum is. Since I is a constant, we can choose it 
to be a vector Le z of fixed length pointing along the 
space fixed unit vector e 2 along the z-axis. When 
n = (sin 9 cos <j), sin 9 sin 0, cos 9) is expressed in polar co- 
ordinates relative to the z axis, Eq. (|A1(I reduces to the 
equations 



dt 
(W_ 
dt 



L 

T 



1 AX 



L sin 29. 



(A2) 
(A3) 



The equation for 9 is easily solved subject to the bound- 
ary condition that 9(t = 0) = 9 : 



tan9(t) = tan^oexp 



1 AI 



4 7 In I a 



(A4) 



Thus, if 2j| > 21, 9(t) -> if < 6> < tt/2 or 0(i) -> tt 
if7r/2<#o<7''asi^oo. This means that n will align 
or anti-align with the angular momentum direction and 
that the angular momentum comes entirely from spin- 
ning parallel to the anisotropy axis with kinetic energy 
L 2 /(2/||). If I\\ < 2_L, 9{t) -> tt/2 for < 9 < tt. In 
this case, n lies in the xy-plane and rotates according 
to Eq. HA2jl. and the kinetic energy is L 2 /(2I±). Thus, 
when angular momentum is conserved but energy is not, 
the rigid body will evolve toward the state with the low- 
est kinetic energy consistent with the constraint of fixed 
angular momentum. 



Here we shortly demonstrate that the antisymmetric 
part of the stress tensor is equivalent to a torque acting 
on the intrinsic angular momentum. The total angular 
momentum of a system with volume V is given by 



L = / (px x v + l)d 3 x 
Jv 



(Bl) 



In case of zero body forces and torques (which might orig- 
inate from applied magnetic and electric fields or grav- 
itation), only surface forces can change the total angu- 
lar momentum. Per definition they are described by the 
stress tensor cr so that 



^ = / x x erdf 
dt 1 



(B2) 

where dV means surface of V . Applying Gauss's theorem 
to the right-hand side results in 



— = / (x x diver + r)d 3 x 

dt Jgy 



where we introduced the torque 
With 



dL 



dv dl N 



, , px x 

dt Jy V at at J 

and the momentum balance in differential form, 
dv 

p- = d 1V< r , 



dl _ 

~dt ~ T ' 



we obtain from Eq. (jB 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



where r is the torque acting on I. 
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